Limits on the WIMP-Nucleon Coupling 
Coefficients from Dark Matter Search 
Experiment with NaF Bolometer 

A. Takeda a '*'\ M. Minowa a ' b , K. Miuchi d , H. Sekiya a , 
Y. Shimizu a , Y. Inoue b ' c , W. Ootani c , Y. Ootuka 6 

^Department of Physics, School of Science, University of Tokyo, 7-3-1, Hongo, 
Bunkyo-ku, Tokyo 113-0033, Japan 

h Research Center for the Early Universe(RESCEU), School of Science, University 
of Tokyo, 7-3-1, Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 

c International Center for Elementary Particle Physics(ICEPP), University of 
Tokyo, 7-3-1, Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 

d Cosmic-Ray Group, Department of Physics, Faculty of Science, Kyoto 
University, Kitashirakawa, Sakyo-ku, Kyoto, 606-8502, Japan 

G Institute of Physics, University of Tsukuba, 1-1-1 Ten'nodai, Tsukuba, Ibaraki 

305-8571, Japan 



Abstract 

We have performed the underground dark matter search experiment with a sodium 
fluoride (NaF) bolometer array from 2002 through 2003 at Kamioka Observatory 
(2700 m.w.e.). The bolometer array consists of eight NaF absorbers with a total 
mass of 176 g, and sensitive NTD germanium thermistors glued to each of them. 
This experiment aims for the direct detection of weakly interacting massive particles 
(WIMPs) via spin-dependent interaction. With an exposure of 3.38 kg days, we de- 
rived the limits on the WIMP-nucleon coupling coefficients, a p and a n . These limits 
confirmed and tightened those derived from our previous results with the lithium 
fluoride (LiF) bolometer. Our results excluded the parameter space complementary 
to the results obtained by Nal detectors of UKDMC experiment. 
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1 Introduction 



Several lines of evidence indicate that the universe contains a large amount of 
nonbaryonic dark matter. The most plausible candidates for this nonbaryonic 
dark matter are weakly interacting massive particles (WIMPs) provided by the 
supersymmetry in the form of the lightest supersymmetric particles (LSPs). 
LSPs can be directly detected through their elastic scattering off target nu- 
clei in a detector [1]. Many experiments have been performed for the direct 
detection of WIMPs. Among them, the DAMA experiment has reported an 
annual modulation signal from the four years' measurement with nine 10 kg 
Nal(Tl) detectors [2]. The parameters reported by the DAMA experiment are 
M WIMP = (52^°) GeV for the WIMP mass and a^ IMP _ p = (7.2±g;f) x 10~ 6 pb 
for the WIMP-proton spin-independently (SI) coupled cross section. On the 
other hand, the CDMS [3] and EDELWEISS [4] experiments have excluded 
most of the parameter space claimed by the DAMA experiment under the as- 
sumptions of standard WIMP interactions and a halo model. Since the cross 
section for the WIMP-nucleus elastic scattering is given by the sum of SI and 
spin-dependent (SD) terms, the dark matter search experiments should be 
performed both via SI and SD interactions. In an SI and SD mixed frame- 
work [5], (Xwimp-p could be as hi g h as °- 6 P b for 3 °- 80 GeV WIMPs in the 
DAMA allowed region under the assumptions of certain parameter sets. Since 
the limits for SD WIMP-proton cross section, <t^ mp obtained by the con- 
ventional method [6] are highly WIMP model dependent, both OyJiMP-p and 
SD WIMP-neutron cross section, o"wJ MP _ n , limits or the limits in terms of the 
fundamental WIMP-nucleon coupling coefficients a p and a n should be derived 

[7]- 

We, the Tokyo group, reported our previous results with the lithium fluoride 
(LiF) bolometer at Kamioka Observatory (2700 m.w.e.) [8], and set the new 
limits in the a p -a n plane. We compared our results with those of the UKDMC 
experiments [10,11], and tightened the UKDMC's limits that are close to the 
DAMA's SD allowed region in Ref. [5]. 

Recently, we developed a sodium fluoride (NaF) bolometer array, and per- 
formed a dark matter search measurement from 2002 through 2003. Since 
23 Na has a larger expectation value of the neutron spins than 7 Li as shown in 
Table 3, we were able to tighten the limits from the previous experiment [8] 
in the a p -a n plane. 

In this paper, we present the results from a dark matter search experiment 
performed at Kamioka Observatory using the NaF bolometer array. 
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2 Experimental set-up 



The present data were collected at the Kamioka Observatory under a 2700 
m.w.e. rock overburden located at Mozumi Mine of the Kamioka Mining and 
Smelting Co. in Kamioka-cho, Gifu, Japan. Details of the underground labo- 
ratory and the experimental set-up are given in Ref. [8]. 

The detector is an array of eight NaF bolometers with a total mass of 176 g. 
Each bolometer consists of 2 x 2 x 2 cm 3 NaF crystal and a neutron transmu- 
tation doped (NTD) germanium thermistor glued to it. Each NTD thermistor 
is biased through a 300 MQ load resistor, and the voltage change across the 
thermistor is fed to the voltage sensitive amplifier with a gain of 10 placed at 
the 4 K cold stage. The signal from the cold stage amplifier is in turn amplified 
by the voltage amplifier with a gain of 4.86 x 10 2 placed at the room tempera- 
ture, and is digitized with a sampling rate of 1 kHz. All the data are recorded 
continuously without any online trigger for a complete offline analysis. 

We refer to each bolometer as from Dl to D8. Five of the NaF crystals are pur- 
chased from OHYO KOKEN KOGYO Co. Ltd., two of them are from RARE 
METALLIC Co. Ltd., and the other is from MATERIAL TECHNOLOGY 
Co. Ltd. We used the NTD thermistors, which are referred to as NTD3 in 
Ref. [12], for the eight bolometers. The concentration of radioactive contam- 
inations in the NaF crystals was previously checked with a low-background 
germanium spectrometer placed at Kamioka Observatory (2700 m.w.e.), and 
was found to be less than 0.1 ppb for U, 1 ppb for Th, and 2 ppm for K. In 
order to eliminate remaining grains of the polishing powder, surfaces of the 
NaF crystals were etched with ultra-pure water in a clean bench of class 100. 

The bolometer array is encapsulated in an inner lead shield with a thickness 
of 2 cm. In the previous measurement, this inner shield was made of over 200 
years old lead. We used a new inner shield with a low-radioactivity lead older 
than 400 years for this experiment. 



3 Measurement and analysis 

From 23rd December 2002 through 24th January 2003, a dark matter search 
measurement was performed. In this work, D6 and D7 showed poor signal to 
noise ratios due to the incomplete thermal conduction between the crystal and 
the thermistor and due to the high electric noises, respectively. We therefore 
analyzed the data obtained with the other six bolometers. 

We performed numerical pulse shaping to all the data using the optimal filter, 
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or the Wiener filter [13]. The energy of each event was determined by the peak 
height of a pulse after shaping. The shaped waveform U (cu) in the Fourier space 
is given by 



u{uj) ~ "W (1) 

where S(u>), C{uj) and N{uj) are Fourier transformations of the ideal signal 
waveform, the measured waveform and the noise, respectively, and the aster- 
isk denotes complex conjugation. Averaged signal of the gamma-ray source of 
137 Cs was used as the ideal signal waveform, while the noises were obtained 
by gathering data while no pulse, whose height above a sufficiently low level, 
exists. A good lineality of this shaping method was confirmed using simulated 
pulses which were obtained by adding the ideal signals of known energy to 
the raw data. Since the signals shaped from the real events have the symmet- 
rical waveform, the asymmetrical waveform from any electric noise events or 
any events that were thought to have occurred within the thermistors were 
eliminated by the waveform analysis. The event selection efficiency was esti- 
mated by analyzing the simulated pulses in the same manner as the real data 
because the deficiency is mainly due to the threshold smeared by the baseline 
fluctuation caused by the noises. Energy resolution of each detector was also 
estimated by a similar method, and listed in Table 1. Energy threshold of each 
detector was set at the energy in which the efficiency was equal to 0.8, and the 
events above this threshold were selected as real bolometer signals. Finally, 
those events which have associated signals in other bolometers within ±10 ms 
were discarded because the expected rate of the coincident WIMP events is 
extremely low. 

Calibrations were carried out using the 60 Co and 137 Cs gamma-ray sources by 
the same manner as described in Ref. [8]. Obtained gains for six bolometers 
are listed in Table 1. Gains and linearities were checked every four or five days, 
and the gains were found to be stable within ±5%. 

Fig. 1 shows the spectra obtained with the six bolometers together with one of 
the spectra (D6) obtained in the previous measurement with the LiF bolometer 
[8]. Live times of the six bolometers are listed in Table 1. The total exposure 
is 3.38 kg days. From this measurement with the NaF bolometer, we obtained 
the same background levels as the previous measurement. 

From the obtained spectra, we derived the limits on the SD WIMP-proton 
elastic scattering cross section, o"wi M p_ p , in the same manner as used in Refs. 
[6,8]. The limits are shown in Fig. 2. The best limit is 27 pb for Mwimp = 
20 GeV. Here we used the astrophysical and nuclear parameters listed in Table 
2. We used the A 2 J( J + 1) values calculated by assuming the odd group shell 
model [18], where A is Lande factor and J is the total spin of the nucleus. 
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Detector 


Dl 


D2 


D3 


D4 


D5 


D8 


Bias Current [nA] 


0.05 


0.07 


0.21 


0.04 


0.04 


0.04 


Gain [mV/keV] 


1.4 


1.8 


0.88 


0.80 


1.3 


0.95 


Threshold [keV] 


14.5 


17.2 


11.6 


37.1 


17.8 


28.0 


Event rate [l(T 3 Hz] 


1.1 


0.51 


0.65 


1.0 


0.95 


1.3 


Resolution (FWHM) [keV] 


1.9 


3.2 


1.4 


6.4 


2.8 


5.4 


Mass [g] 


21.5 


21.5 


21.7 


21.6 


21.5 


23.4 


Live Time [days] 


26.0 


25.2 


26.0 


26.0 


25.5 


25.9 



Table 1 

Detector responses of the six bolometers. 



Dark matter density 


PD = 


0.3GeVc- 2 /cm 3 


Velocity distribution 


Maxwellian 


Velocity dispersion 


v = 


220 kms" 1 


Escape velocity 


V esc = 


= 650kms~ 1 


Earth velocity 


V E = 


217 kins" 1 



Spin factor ( 23 Na) A 2 J( J + 1) = 0.041 
Spin factor ( 19 F) A 2 J( J + 1) = 0.647 



Table 2 

Astrophysical and nuclear parameters used to derive the exclusion limit. 

As discussed previously, these cr^ MP limits are fraught with potentially 
significant WIMP model dependence. We then derived the model-independent 
limits on 0v^ MP _ p and 0wiMP-n m the same procedure as described in Ref. 
[7]. The SD WIMP-nucleus cross section, cr^ MP _ N , is given by the following 
equation: 



°"WIMP-N - 4G F /i WIMp _ N X (dp (S'p(N) ) + On^OM))) 2 J~ , (2) 

= 4Gp/X 2 VIMp _ N X (-y/Cp(N) ± \J Cn(N)) 5 (3) 

where Gf is the Fermi coupling constant, /xwimp-n is the WIMP-nucleus re- 
duced mass, a p and a n are the WIMP-proton and WIMP-neutron coupling 
coefficients, respectively, (S p ^) and (S'n(N)) are the expectation values of the 
proton and neutron spins within the nucleus N, respectively, Cp(N) and 
are the proton and neutron contributions to the total enhancement factor 
of nucleus N, respectively. The sign in the square of Eq. (3) is determined 
by the relative sign of a p (S p ^) and a n (S n ^). (S p ^) and (S n ^) calculated 
with recent shell models are shown in Table 3. From the WIMP-nucleus cross 



5 



section limit, ct^mp_ n , which is directly set by the experiments, the respec- 
tive limits on WIMP-proton and WIMP-neutron cross section, ^^mp„ p (n) or 
^wiMP-nfN)) are given by assuming that all events are due to WIMP-proton 
and WIMP-neutron elastic scatterings in the nucleus, such as: 



SDlim _ SDlim PWIMP-p / ^p(N) . / 

°WIMP-p(N) — ^WIMP-N"! I "T^sd" j ' V I 



..2 /r~<SD 
PWIMP-p / u p ' 

A^WIMP-N \ ^p D 



-1 



2 /r SD \ _1 

SDlim _ SDlim ^WIMP-n f °n(N) \ /_>. 

°WIMP-n(N) _ ^WIMP-N - 2 \ q SD J j \°) 



A*WIMP-N 



where /xwimp-p an d /-twiMP-n are the WIMP-proton and WIMP-neutron re- 
duced mass, respectively, and and C„ D are the enhancement factors of 
proton and neutron themselves, respectively. Ratios of the enhancement fac- 
tors calculated by the shell models are shown in Table 3. 

The obtained model-independent limits on o"^ M p_ and 0v^ MP _ n are shown 
in Fig. 3. We then derived the limits in the a p -a n plane from them. According 
to Ref. [7], the allowed region in the a p -a n plane is defined as: 



± ; ^ 1 , (6) 



/_.SD lim / _SD lim I ~~ 24Gl „t> 

.V^WIMP-pCN) V ^WIMP-n(N)/ ^^F^WIMP- 



where the relative sign inside the square is determined by the sign of (^(n))/ (<Sn(N)), 
and the small proton-neutron mass difference is ignored. Obtained limits in 
the a p -a n plane for WIMPs with mass of 10 GeV, 50 GeV and 100 GeV are 
shown in Fig. 4. The limits on the respective coupling coefficients are, for 
instance, \a p \ < 16 and \a n \ < 68 for M W imp = 50 GeV. 



In Fig. 4, we compared our results with those of our previous results [8] and 
those of UKDMC experiment [11]. Though our 0wiMP-p limits derived from 
this measurement are almost the same as those derived from the previous 



Isotopes 


Abundance (%) 


J 


('S'p(N)) 


(£n(N)) 


°p(N)/°P 


U n(N)/ n 


7 Li 


92.5 


3/2 


0.497 


0.004 


5.49 x 10- 1 


3.56 x 10- 5 


19p 


100 


1/2 


0.441 


-0.109 


7.78 x 10- 1 


4.75 x 10~ 2 


23 Na 


100 


3/2 


0.248 


0.020 


1.37 x 10- 1 


8.89 x 10~ 4 


127j 


100 


5/2 


0.309 


0.075 


1.78 x 1Q- 1 


1.05 x 10~ 2 



Table 3 

Values of (S* p ( N )), (S n ^), C^^/C^, C^^/C^ for various isotopes. Values for 7 Li 
and 19 F are taken from Ref. [19] and those for 23 Na and 127 I are taken from Ref. [20] 



6 



results with the LiF bolometer, we improved the limits in the a p -a n plane 
significantly. This is because 23 Na has a larger expectation value of the neutron 
spins than 7 Li. It is clearly seen that our results obtained with LiF or NaF 
detectors set limits complementary to the results obtained with Nal detectors 
of UKDMC experiment because the spin factor of 19 F are large and the sign of 
(-S'p(N))/('S'n(N)) for 19 F is opposite to those for 7 Li, 23 Na, and 127 I. As described 
in Ref. [8], we did not compare our results with the SD annual modulation 
allowed regions shown in Fig. 2 [5] because the DAMA group had not shown 
the results from 23 Na and 127 I independently. 

The remaining background in the energy region between 10 and 90 keV re- 
stricting the current sensitivity is supposed to be due to the uranium and 
thorium contamination in a crystal holder made of the OFHC copper. A use 
of the crystal holder with a sufficiently radio-pure material such as polyte- 
trafluoroethylene (PTFE) and active shields between the crystals and their 
holder would help to reduce the background from the holder and the other 
outside materials. 



4 Conclusion 

We have performed the dark matter search experiment at Kamioka observa- 
tory with the NaF bolometer. With a total exposure of 3.38 kg days, we have 
confirmed and tightened the limits in the a p -a n plane obtained by the previous 
results with the LiF bolometer. For the WIMPs with a mass of 50 GeV, we 
have successfully excluded more than half of the parameter region in the a p -a n 
plane allowed by the previous results with the LiF bolometer. Our results set 
limits complementary to the results obtained by Nal detectors of UKDMC 
experiment. Since the limits from the UKDMC experiments are close to the 
DAMA's SD allowed region in Ref. [5], we should also be able to compare our 
limits in the a p -a n plane if DAMA publishes relevant data to calculate these 
parameters. It is worth noting that this result has been obtained by a single 
detector with both 23 Na, which has relatively large C^^/C^ , and 19 F, which 
has the largest C p (N)/C p D , in the detector material. 
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Fig. 1. Low energy spectra obtained with the six NaF bolometers are shown. One 
of the spectra (D6) obtained in the previous measurement with the LiF bolometer 
[8] is also shown on top. 




Fig. 2. NaF combined 90% C.L. 0\yiMP-p ^ m ^ s as a function of Mwimp are shown 
in a thick solid line. Limits from our previous results with the LiF bolometer at 
Kamioka Observatory [8] and from other experiments [11,14,9,15,16,17] are shown 
in a thick dash-dotted and thin dotted lines, respectively. DAMA's allowed region 
investigated in a mixed coupling framework shown as "case c" in Ref. [5] is shown 
in a closed contour. 




10° 10 1 10 2 10 3 10° 10 1 10 2 10 3 
Wwimp [GeV] M W imp [GeV] 

Fig. 3. Model independent 90% C.L. limits on spin-dependent WIMP-proton (a) 
and WIMP-neutron (b) cross sections as a function of Mwimp- The limits have 
been derived by the same manner as described in Ref. [7]. The limits from Na and 
F are shown in thin dotted lines. The combined limits are also shown in thick solid 
lines. 




-40 -20 20 40 

a P 

Fig. 4. Limits in the a p -a n plane: (a) For Mwimp = 10 GeV, (b) for 
-^wimp = 50 GeV, and (c) for Mwimp = 100 GeV. Inside of thick solid lines are 
allowed region by this work. The limits from our previous results with the LiF 
bolometer at Kamioka Observatory and UKDMC data[ll] are also shown for com- 
parison in thick dotted lines and thin dashed lines, respectively. 



